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Abstract 

The Airborne Cloud-Aerosol Transport System (ACATS) is a multi-channel Doppler 
lidar system recently developed at NASA Goddard Space Flight Center (GSFC). A 
unique aspect of the multi-channel Doppler lidar concept such as ACATS is that it is also, 
by its very nature, a high spectral resolution lidar (HSRF). Both the particulate and 
molecular scattered signal can be directly and unambiguously measured, allowing for 
direct retrievals of particulate extinction. ACATS is therefore capable of simultaneously 
resolving the backscatter/extinction properties and motion of a particle from a high 
altitude aircraft. ACATS has flown on the NASA ER-2 during test flights over California 
in June 2012 and science flights during the Wallops Airborne Vegetation Experiment 
(WAVE) in September 2012. This paper provides an overview of the ACATS method 
and instrument design, describes the ACATS retrieval algorithms for cloud and aerosol 
properties, and demonstrates the data products that will be derived from the ACATS data 
using initial results from the WAVE project. The HSRL retrieval algorithms developed 
for ACATS have direct application to future spacebome missions such as the Cloud- 
Aerosol Transport System (CATS) to be installed on the International Space Station 
(ISS). Furthermore, the direct extinction and particle wind velocity retrieved from the 
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ACATS data can be used for science applications such as dust or smoke transport and 
convective outflow in anvil cirrus clouds. 

1.0 Introduction 

Current uncertainties in the role of aerosols and clouds limit our ability to accurately 
model the Earth’s climate system and to predict climate change. There are several 
different types of lidar systems that can be used to measure cloud and aerosol properties 
and motion. Cloud-aerosol lidars measure the elastic backscatter from molecules and 
atmospheric particulates to resolve vertical profiles of spatial and optical properties of 
clouds and aerosols. The two most common elastic backscatter lidar techniques are 
standard backscatter lidars and high spectral resolution lidars (HSRL). The data provided 
by these lidar systems are essential to investigations of cloud and aerosol properties for 
numerous reasons. The vertical structure of cloud and aerosol layers resolved by lidar 
systems cannot be accurately obtained from passive satellite or passive airborne sensors. 
Furthermore, thin cloud optical depths are often below the detection limits of millimeter 
cloud radar systems (Comstock et al. 2002). In situ instruments can provide critical 
measurements of cloud and aerosol microphysical properties. However, they do not 
easily provide vertical profiles of these measurements and can alter the physical 
properties of the particles (Jensen et al. 2009; Zhao et al. 2011). Information obtained 
from cloud-aerosol lidar systems can improve knowledge of cloud and aerosol properties, 
which in turn advance parameterizations and reduce the uncertainties introduced in 
GCMs. 

Standard elastic backscatter lidars are the least complex and most common lidar 
systems used to study vertical profiles of cloud and aerosol properties. Ground-based and 
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airborne systems have been used in numerous field campaigns over the past few decades. 
In the last decade, as laser transmitters have become more reliable, the first space-based 
elastic backscatter lidar systems were designed and launched. The Geoscience Laser 
Altimeter System (GLAS; Spinhirne et al. 2005) was launched in January 2003 and the 
Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations project (CALIPSO; 
Winker et al. 2009) was launched in April 2006. These lidar systems fundamentally 
measure vertical profiles of attenuated total backscatter, without separation of particulate 
(Mie) and molecular (Rayleigh) scattering. There have been many methods developed to 
retrieve the particulate extinction and particulate backscatter coefficients from a cloud- 
aerosol lidar return signal. One technique is an inversion using standard backscatter lidar 
data developed by Femald et al. (1972) and Klett (1981; 1985). The Klett or Femald 
method makes it possible to solve the standard lidar equation by assuming a ratio of 
aerosol extinction to aerosol backscatter coefficients, referred to as the lidar ratio, is 
known and constant throughout a particulate layer. This assumption reduces the number 
of unknowns in the system to one. This method is commonly used to retrieve particulate 
extinction and backscatter coefficients from standard backscatter lidars such as CALIPSO 
(Young and Vaughan 2009) and the Cloud Physics Lidar (McGill et al. 2002). The lidar 
ratio (units of sr) is highly dependent on the optical and microphysical properties of 
atmospheric layer being measured. The lidar ratio typically varies from about 10 to 50 sr 
for tropospheric clouds (Del Guasta et al. 2001; Seifert et al. 2007; Yorks et al. 2011a) 
and from about 20 to 80 sr for aerosol particles (Ackermann 1998). For cloud and 
aerosol layers with an optical depth greater than 0.30, a 30 percent error in the assumed 
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lidar ratio can lead to an error in the extinction retrieval from elastic backscatter lidar 


systems greater than 50 percent (Young et al. 2014). 

Another method for retrieving the particulate backscatter and extinction coefficients 
from a lidar signal is a HSRL, which is based on the use of two measured profiles instead 
of only one. The HSRL technique utilizes the difference in spectral distribution of the 
molecular and particulate backscattered signals (Fiocco et al. 1971; Shipley et al. 1983; 
Grund et al. 1991). High spectral resolution optical filters are required to separate the 
particulate contribution from the molecular backscatter and resolve particulate extinction 
and backscatter coefficients independently with no assumption about the lidar ratio 
required. Only a few HSRL instruments have been successfully developed and operated 
to measure cloud and aerosol optical properties from ground or aircraft platforms. The 
iodine filter method (Piironen and Eloranta 1994) is the preferred method for HSRL 
systems to date. Recently airborne HSRL systems that employ iodine filters have been 
implemented and demonstrated on the NASA King Air (B-200) research aircraft (Hair et 
al. 2008) and the German Aerospace Center (DLR) Falcon research aircraft (Esselborn et 
al. 2008). However, a caveat of the iodine filter technique is that the actual particulate 
backscatter spectral broadening is not measured but inferred from the total and molecular 
backscatter. The backscattered signal also contains additional information that is 
imparted in the scattering process, such as the Doppler shift caused by the mean velocity 
of the particulate. 

Doppler wind lidars use the frequency shift imparted on atmospheric aerosols and 
molecules to determine vertical profiles of the horizontal wind speed and direction. 
Providing these measurements on a global scale can progress understanding of 
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atmospheric dynamics and improve numerical weather predictions (Baker et al. 1995). 
The two most common types of pulsed Doppler wind lidar systems are coherent 
(heteodyne) detection and direct (incoherent) detection. Coherent Doppler lidars use a 
heterodyning technique that mixes a pulsed lidar signal with a second laser signal to 
produce a beat frequency that is related to the Doppler shift. The second continuous laser 
beam is usually a local oscillator offset in frequency (Hall et al. 1984; Huffaker et al. 
1984). Direct-detection lidars directly measure the frequency shift of the return signal 
using a high spectral resolution filter, such as a Fabry-Perot interferometer or etalon, and 
operate at shorter wavelengths than coherent systems (Benedetti-Michelangeli et al. 
1972; Chanin et al. 1989; Gamier and Chanin 1992; Gentry and Korb 1994). One direct- 
detection method, termed multichannel (MC) by McGill and Spinhime (1998), measures 
the Doppler shift by imaging the etalon fringe pattern onto a multiple element detector 
(Abreu et al. 1992; Fischer et al. 1995) The MC direct-detection concept requires the 
etalon transmission function to be aligned with the laser wavelength. This method was 
demonstrated by McGill et al. 1997a for a ground-based lidar developed at the University 
of Michigan. 

A MC direct-detection Doppler lidar system capable of resolving the Doppler shifts 
inherent to atmospheric motions can simultaneously provide information about both the 
scattering intensity and the motion of the particle. Such an instrument was recently 
developed at NASA Goddard Space Flight Center (GSFC) called the Airborne Cloud- 
Aerosol Transport System (ACATS). ACATS is the first lidar system to simultaneously 
measure cloud/aerosol properties and wind from an airborne platfonn. The instrument 
has flown on the NASA ER-2 during test flights over California in June 2012 and as part 
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of the Wallops Airborne Vegetation Experiment (WAVE) in September 2012. A 
description of the ACATS instrument design is provided, which includes details of the 
optical and mechanical components of the subsystems as well as the software that 
autonomously controls the instrument operation. This work advances the effort of McGill 
et al. 1997a and McGill et al. 1997b by demonstrating the retrieval algorithms for HSRL 
direct measurements of cloud and aerosol optical properties (i.e. extinction) that can be 
applied to future space-based HSRL missions. This study also presents initial ACATS 
HSRL results and data products from the WAVE campaign. 

2.0 ACATS Method and Instrument Description 

2.1 ACATS Methodology 

The ACATS instrument is a multi-channel (MC) Doppler lidar system built for use on 
the NASA ER-2 high altitude aircraft. The MC technique passes the returned 
atmospheric backscatter through a single etalon and divides the transmitted signal into 
several channels (wavelength intervals), which are measured simultaneously and 
independently (Figure 1). The resulting aerosol spectral distribution is then compared to 
the outgoing laser distribution to infer the Doppler shift, as demonstrated in Figure 2a. 
Subsequent measurements of the atmospheric scattered light will reveal a wavelength 
offset that is proportional to the Doppler shift and directly related to the velocity of the 
scattering particles (Figure 2b). The basic concept is summarized in Figures 1 and 2. 
The MC method was demonstrated using the ground-based University of Michigan 
Doppler lidar (McGill et al. 1997a; McGill et al. 1997b). 

A unique aspect of the MC Doppler lidar concept such as ACATS is that it is also a 
HSRL. Both the particulate and molecular scattered signal can be directly and 
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unambiguously measured since the broad Rayleigh-scattered spectrum is imaged as a 
nearly flat background, illustrated in Figure 2c. The integral of the aerosol-scattered 
spectrum is analogous to the aerosol measurement from the typical absorption filter 
HSRL technique, providing exactly the same pieces of information as a standard FISRL 
(Figure 2d). While previous ground-based MC systems have been built and operated 
(Benedetti-Michelangeli et al. 1972; Abreau et al. 1992; McGill et al. 1997a), there has 
been no airborne demonstration of the technique and the method has not been used to 
derive HSRL cloud and aerosol properties. 

2,2 ACATS Instrument Description 

The ACATS instrument is composed of three main subsystems; laser transmitter, 
telescope, and receiver optics. A picture of the ACATS instrument fully assembled, with 
the receiver and telescope subsystems, is shown in Figure 3. A list of the ACATS 
instrument parameters is provided in Table 1. The instrument also includes a 
heating/cooling loop to provide stable thermal operation of the laser. 

The frequency characteristics of pulsed lasers have recently been advanced due to the 
development of direct detection Doppler lidars and HSRLs. These techniques impose 
further requirements compared to standard backscatter lidars, such as lasers that are 
single frequency on a single pulse basis and more stable in time (central frequency drift 
of less than 1 MHz per minute). An injection-seeded, pulsed Nd:YAG laser was 
developed for the TWiLiTE instrument (Hovis et al. 2004) that achieves these frequency 
characteristics. This laser was later replicated for the ACATS instrument and provides a 
narrow wavelength distribution suitable for resolving the small frequency shifts due to 
the Doppler effect. The laser operates at an output power of about 10 mJ per pulse and 
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repetition rate of 250 Hz at 532 nm and is designed for use in the low-pressure 
environment of high-altitude aircraft. 

The ACATS telescope employs a rotating holographic optic element (HOE) to fit the 
small volume envelope of the ER-2 superpod and to enable vector wind measurements, 
which requires more than one viewing direction (Figure 3 c). The telescope system is set 
for 45 degree off-nadir viewing and rotates on a bearing to permit step-stare operation. 
The number of scan angles (up to 8) and dwell time at each scan angle is controlled by 
software and can be modified before flight. A schematic of the optical design is presented 
in Figure 4. As the telescope rotates, the optical alignment changes and may lead to a 
loss in return signal if not corrected. A procedure that steps the telescope position using 
piezoelectric actuators and scans for the largest return signal is run during flight to 
determine the optical alignment at each scan position. The 8-inch diameter telescope is 
also fiber-coupled to the receiver subsystem to provide greatest flexibility. 

The primary difference between a lidar system capable of only measuring total 
backscatter intensity (e.g., CALIOP or CPL) and an instrument that directly measures the 
particulate extinction and Doppler shift, such as ACATS, lies in the receiver subsystem 
(Figures 3b; 4). The heart of the ACATS receiver system is an etalon that provides the 
spectral resolution needed for the HSRF measurement and also to resolve the Doppler 
shift inherent in the backscattered signal. Backscattered light collected by the telescope 
is passed through the etalon and an image of the etalon fringe pattern is created. A 
bandpass filter is used in tandem with the etalon to reject background sunlight, permitting 
daytime operation. The optical gap of the etalon is 10 cm with an operational diameter of 
35 nun and plate reflectivity of 85%. As with any MC system, it is critical to maintain 
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the symmetry and shape of the etalon fringe pattern to avoid uncertainty in the 
measurement. A digital etalon controller was developed by Michigan Aerospace 
Corporation in which piezoelectric actuators control the etalon electronics to position and 
maintain the plate parallelism. Considerable work was performed to create autonomous 
flight software that maintains the etalon alignment over the entirety of an ER-2 flight. 
The signal transmitted by the etalon is then passed to the detector subsystem. 

A holographic circle-to-point converter optic (McGill et al. 1997c; McGill and 
Rallison 2001) is placed in the focal plane to provide the spectral detection. The circle- 
to-point converter simplifies hardware requirements, improves efficiency of measuring 
the spectral content in the fringe pattern, and allows ACATS to utilize photon-counting 
detection. The holographic optic is coupled to a Hamamatsu H7260 linear array detector, 
which utilizes back-end electronics developed by Sigma Space Corporation to permit 
photon-counting detection at count rates in excess of 50 MHz. The ACATS receiver 
images ~1.2 orders over 24 detector channels. The ACATS etalon parameters result in a 
measurement dynamic range of -400 m/s, more than sufficient for typical atmospheric 
motions. 

An autonomous multi-channel data system is the final component of the instrument 
and was based entirely on work completed by Sigma Space Corporation in support of the 
CPL, UAV-CPL, and TWiLiTE lidars. The basis for the data system, the Advanced 
Multichannel Scaler (AMCS) card, was first applied in the ER-2 CPL instrument. The 
data acquisition software is included in the data system and has its heritage in the CPL 
and UAV-CPL instruments. An important aspect of the ACATS data system, as 
developed for CPL and UAV-CPL, is the ability to downlink data in real-time from the 
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aircraft using the onboard air and navigation payload server. The data system also 
incorporates a Novatel model OEMV-3RT2i GPS receiver and OEM-IMU-H58 inertial 
unit to enable accurate correction for platform motion. The Novatel system provides 
greater than 20 Hz update rates with 2 cm/s velocity accuracy. The raw ACATS data file 
consists of photon counts at each horizontal record (1 sec), range bin (30 m) and detector 
channel, which is then converted to atmospheric parameters such as backscatter and 
extinction coefficients. 

2.3 ACATS Calibration Procedures 

Several calibration parameters are required to accurately retrieve the wind velocity, 
aerosol and molecular backscatter from the ACATS data. These include normalization 
constants, instrument defect function parameter, and detector nonlinearity. The 
illumination and sensitivity of each detector channel are not the same, necessitating 
normalization constants to compensate. The detector normalization coefficients are 
detennined using a white-light source to illuminate the telescope while the receiving 
optics remains unchanged. These nonnalization constants describe the relative response 
of the detector to broad bandwidth illumination. 

The alignment of the circle-to-point converter (HOE) and Fabry-Perot fringe pattern 
also must be characterized. Each ring in the circle-to-point converter represents a 
detector channel. Since the circle-to-point converter and etalon are manufactured 
separately, a ring can have a dissimilar centricity and diameter compared to the fringe 
pattern projected onto it, resulting in signal loss to the corresponding detector channel. 
To complicate matters, this loss of signal can vary in each channel. In the case of 
ACATS the outer rings (higher detector channels) of the circle-to-point converter are not 
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perfectly concentric with the fringe pattern, requiring normalization constants to 
compensate. The nonnalization coefficients are determined using the peak transmission 
of the etalon calibration data in each channel. Assuming perfect alignment in all 
channels, the peak transmission will remain constant as the signal is stepped through all 
detector channels. Thus, the ACATS channel with the highest transmission represents 
the best alignment, allowing all other channels to be normalized to the “best aligned” 
channel. These normalization constants describe the relative signal loss of the detector 
channel due to alignment imperfections. 

To characterize the instrument defect parameter, an etalon calibration procedure has 
been developed for ACATS similar to the one outlined in McGill et al. (1997a). The 
etalon transmission equation as a function of detector channel (j) is expressed as (McGill 
1996): 


T(MJ) 


XA, 

( 1=0 
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2 Jin 


AA j 

• + • J 


\ ^A-FSR 
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sine 


FSR / 
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FSR / 


(2.3.1) 


where AXfsr is the free spectral range and is defined as the change in wavelength 
necessary to change the order of interference by one. The free spectral range can also be 
represented by the number of channels necessary to change the order of interference by 
one, Nfsr- The function A n is defined as: 


A„ =21- 


£ Wl -R\ 


n -4 jt 2 n 2 A d ^ Aq 


R n e 


(2.3.2) 


1 -R) U + RJ 

where £ is the loss of light due to absorption or scattering by the etalon plates and R is 
the plate reflectivity. The etalon transmission equation (2.3.1) is for an idealized etalon. 
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A real etalon function will be broadened by several effects, such as plate bowing, 
microscopic plate defects, detector broadening, and off-axis aberrations. 

For the purpose of this study, it is sufficient to use an instrument defect parameter 
(Ad D ) to represent the etalon broadening effects and tune the etalon model so that it 
matches the measured AC ATS spectral response. There are two important assumptions in 
determining the ACATS defect parameter. First, the defect parameter varies with detector 
channel to account for the variability of the etalon finesse with channel. It is also assumed 
that any broadening effects, and thus the etalon defect parameter, will follow a Gaussian 
distribution. The ACATS defect parameter is then determined by a calibration procedure 
similar to the one demonstrated in McGill et al. (1997a). Software runs a calibration 
procedure at least once per flight that varies the etalon gap using piezoelectric actuators. 
Varying the etalon gap moves the interference fringe pattern across the detector in 128 
small steps, sampling nearly 3 orders (42 points per order). One can then determine the 
defect parameter for each channel by performing a least-squares fit to match the modeled 
etalon transmission function to the ACATS measured etalon response function using a 
similar technique to McGill et al. 1997a. The light source used to measure the ACATS 
etalon response is the same laser that is used for atmospheric measurements. 
Additionally, the calibration technique automatically compensates for any uncertainty in 
computing the laser bandwidth, since the laser width follows a Gaussian distribution 
similar to the etalon broadening term. 

The measured ACATS spectrum can become distorted due to detector dead time and 
must be compensated for. All lidar systems that employ photon-counting detection 
experience this effect, which is a limitation on the number of photons that can be counted 
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in a given time interval. For ACATS, the large near- field return pushes the detector into 
a nonlinear counting region. The nonlinear effects for this type of detector can be 
quantified by a detector dead time coefficient. This coefficient represents the fact that 
only one photon event can be counted at once, and the detector system has a certain time 
delta, or dead time, before it can count another. A typical Hamamatsu linear array 
detector, such as the one employed in ACATS, has a discriminator dead time of 65 to 75 
ns for a discriminator maximum count rate on the order of 15 MHz. To improve this 
performance, the ACATS Hamamatsu linear array detector is customized with a 
discriminator built by Sigma Space Corporation under Small Business Innovative 
Research (SBIR) funding that has a shorter discriminator dead time. This permits 
photon-counting detection at count rates in excess of 40 MHz before there is a 10% 
reduction in observed count rate. The ACATS detector rarely experiences count rates 
higher than 10 MHz in atmospheric bins below 17 km (assuming an ER-2 altitude greater 
than 19 km). Therefore, the detector dead time coefficient is less than 1.05 for 99.5% of 
atmospheric bins with the exception of the near-field return. 

3.0 Development of ACATS Retrieval Algorithms 

ACATS provides data products similar to other cloud-aerosol lidars, HSRL systems, 
and Doppler wind lidars. The system is currently set for 45 degree off-nadir viewing and 
the telescope rotates to allow for two orthogonal line-of-sight (LOS) wind measurements, 
which are then used to compute vertical profiles of horizontal wind velocity and direction 
within particulate layers. The ACATS retrieval algorithms and data products for the 
horizontal wind velocity will be presented at a later date. This paper focuses on two 
types of aerosol/cloud products available from ACATS data that are directly applicable to 
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the ISS CATS instrument. Standard backscatter products are computed similar to CPL 
and CALIPSO (McGill et al. 2007). HSRL products are produced at courser resolutions 
(450 m vertical and 5 km horizontal), but include direct retrievals of attenuated 


particulate backscatter, optical depth, as well as particulate extinction and backscatter 
coefficients. These products are similar to those produced by other HSRL systems. 

3.1 Development of Standard Backscatter Algorithms 

If the measured ACATS photon counts are summed over all channels as to neglect the 
spectral information provided by the etalon, vertical profiles of total backscatter can be 
retrieved from ACATS data. Similar to a standard backscatter lidar system (i.e. 
CALIOP), this total signal is composed of both the particulate scattering and molecular 
scattering. The standard lidar equation can be regrouped and solved for the attenuated 
total backscatter (ATB or y), which has units of km 1 sr 1 and is defined as: 


The molecular backscatter coefficient ((3m) is detennined from Rayleigh scattering theory 
(Tenti et al. 1974; Young 1981) and is proportional to atmospheric density. Furthermore, 
the molecular extinction coefficient (om) is resolved from the molecular backscatter 
coefficient though the relationship a M (k) = IKiffiy) * ( 8/3 )xc. The ACATS standard 
ATB is computed using the standard lidar equation and calibrated by normalizing the 
signal to the molecular backscatter profile at high altitudes where aerosol loading is 
weakest (Russell et al. 1979; Del Guasta 1998). This calibration technique is the well- 
accepted method of calibrating backscatter lidar signals and is used in CALIPSO and 
CPL retrievals (McGill et al. 2007). ACATS cloud and aerosol layer boundaries are 
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determined using a similar method to CPL (Yorks et al. 2011b). The advantage of using 
this retrieval scheme is that the particulate layer properties can be obtained at higher 
resolutions, both vertically and horizontally, than using the HSRL retrieval algorithms. 
Therefore, this “standard” lidar method is used to compute ACATS attenuated total 
backscatter, as well as cloud and aerosol layer boundaries at a vertical resolution of 40 m 
and horizontal resolution of 400 m (2 sec). 

3.2 Development of HSRL Algorithms 

The ACATS HSRL retrieval algorithms are unique and different compared to the 
algorithms of current iodine filter HSRL systems (Hair et al. 2008). The inclusion of an 
etalon in the ACATS instrument design results in a more complicated ACATS lidar 
equation compared to the standard lidar equation and iodine filter HSRL equations. The 
etalon transmission function (Equation 2.3.1) is convolved with the standard backscatter 
lidar equation to yield the expression for the number of photon counts detected per 
channel (j), as derived in McGill 1996: 
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(3.2.1) 

The first term represents the instrument parameters and the definitions of individual 
parameters are shown in Table 2. The second term contains the laser broadening (AXl), 
molecular broadening (AX M ), and the atmospheric physics. The attenuated particulate 


backscatter (a) and attenuated molecular backscatter (to) are expressed as: 
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The Doppler shift is characterized by the second part of the third term, where Ulos is the 
LOS wind velocity in ms' 1 . The attenuated particulate backscatter, attenuated molecular 
backscatter, and LOS wind velocity are the three unknown variables in Equation 3.2.1. 
Since there are 24 detector channels, the ACATS system is an over-determined set of 
equations. These three unknowns are determined using a method developed by McGill et 
al. (1997b). First, the ACATS lidar equation (Equation 3.2.1) is linearized by expanding 
the relevant variables in a Taylor series. The equation is then written in matrix fonn: 
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This equation can also be written as: 


A7V = GAx 


(3.2.5) 


An iterative weighted least-squares fitting technique is employed to resolve these three 
parameters and their corresponding uncertainty, in which the solution is: 
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where W is the weighting matrix and G is the generalized matrix to be inverted. The 
solution for the molecular and particulate signals are linear, but non-linear for the 
Doppler shift. This least-squares fit method was tested and proven by McGill et al. 
(1997b) to retrieve the horizontal wind velocity. This work advances the effort of McGill 
et al. (1997a) and McGill et al. (1997b) by developing HSRL retrievals of cloud and 
aerosol properties. The first step is to compute the molecular backscatter coefficient ((3m) 
and two-way transmission (Tm 2 ) from Rayleigh scattering theory and meteorological data 
from a nearby radiosonde. The definition for the attenuated molecular backscatter (Eq. 
3.2.2) can be rewritten in terms of the two-way transmission, corrected for the slant path, 
and solved for the two-way particulate transmission (T p 2 ): 


T;(r) = 




-ICOS0 


(3.2.7) 


Therefore, the two-way particulate transmission can be determined without making 
unnecessary assumptions about the lidar ratio, as in the Klett or Femald method (Femald 
et al. 1972; Klett 1981, 1985). Once T p 2 is known, the definition of the attenuated 
particulate backscatter (Eq. 3.2.3) can be rewritten and used to directly retrieve the 
particulate backscatter coefficient (|3p): 


/3 P (x,r) 


a(jt,r ) 

Z2(r)7?(r) 


(3.2.8) 


The particulate optical depth is then: 

T,(r)=-lln[r, 2 (r)] 


(3.2.9) 


The particulate extinction coefficient (op) is directly retrieved using the equation: 
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(3.2.10) 


cr P (r) = 


dTp(r) 

dr 


and the particulate lidar ratio is: 


S P (r) = 


Qp(r ) 
P P (r) 


(3.2.11) 


This method is used to compute profiles and layer-integrated values of the 
aforementioned variables at a vertical resolution of 450 m and horizontal resolution of 5 
km (25 sec). Their corresponding uncertainties are computed using propagation of errors. 
If high-resolution optical properties are desired, the directly retrieved lidar ratio can be 
utilized as a parameterization to compute high-resolution optical properties using the 
Klett or Fernald method. 

4.0 Initial Results from WA VE Campaign 

During the period of 9 to 27 September 2012, ER-2 aircraft flights were 
conducted out of Wallops Island, VA as part of the WAVE project. These flights were 
planned over land, targeting specific land and vegetation surfaces with a scientific 
objective of simulating Ice, Cloud and land Elevation Satellite 2 (ICESat-2) data using 
the Multiple Altimeter Beam Experimental Lidar (MABEL; McGill et al. 2013). AC ATS 
was a payload on a total of 13 ER-2 flights, which included observations of thin cirrus 
clouds, and smoke layers. During these flights, software directed the ACATS telescope to 
rotate counter-clockwise to four look angle positions denoted by azimuth angle relative to 
the aircraft nose: 0° (fore), 90° (right), 180° (aft), and 270° (left). At each look angle, the 
dwell time was set for 60 seconds. The WAVE campaign represents the first science 
flights for the ACATS instrument in which the telescope rotated and more than one look 
angle was used. Due to limited time before the project, the telescope alignment was 
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optimized only at the 270-degree look angle. The telescope alignment for the other three 
look angles was performed in the field using the new and untested in-flight telescope 
alignment procedure. Portions of flights, and in some cases entire flights, were used to 
test and refine the etalon calibration procedure and telescope alignment. Furthermore, 
only two look angles were used for some flights if proper telescope alignment was not 
achieved at all four look angles. An example of the photon counts summed across all 24 
detector channels at each of the four look angles from the 26 September 2012 flight is 
shown in Figure 5 and demonstrates the ability of ACATS to observe cirrus clouds 
(between 10 and 12 km) at multiple look angles. Overall, ACATS collected science data 
with high signal-to-noise ratio (SNR) in at least one look angle during 8 of the 13 total 
flights. The telescope alignment and LOS wind retrievals will be improved before future 
ACATS flights. This study will focus on ACATS retrievals of cloud and aerosol 
properties from the WAVE project, particularly those at the 270-degree look angle and 
high quality data from the other look angles. 

There were several flights during WAVE in which ACATS collected quality data 
at multiple look angles. Perhaps the best ACATS perfonnance was on the 26 September 
ferry flight back to Palmdale, CA when all four look angles were well aligned. Figure 6 
shows the 532 nm ATB (km' 1 sr' 1 ) computed using the standard method (a), the 
Attenuated Particulate Backscatter (km' 1 sr' 1 ) using the HSRL method (b), and the 
directly-retrieved Particulate Extinction Coefficient (km' 1 ) at the 0 degree look angle (c) 
for the flight on 26 September 2012. Clearly visible in these images are cloud layers 
observed by ACATS as the ER-2 flew over the Ohio River Valley (20:28:05 to 21:30:00 
UTC) and over North Dakota (about 00:24:10 UTC). ACATS also measured a large 
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smoke plume (00:24:10 to 02:10:00 UTC) that extended as high as 6 km over Montana. 
The images in Figure 6 demonstrate the typical ACATS cloud and aerosol data products. 
The extinction and backscatter values are typical for cloud and smoke layers and appear 
to be similar across retrieval methods. 

The ACATS telescope alignment on the 14 Sep. flight at the 270 degree look 
angle was the best for the entire campaign, making it a good case to assess biases in the 
two retrieval methods. Figure 7 shows the 532 nm ATB computed using the standard 
method (a) and using the HSRL method (b). The latter is essentially a + oo. Cirrus 
clouds between 9 and 13 km are observed throughout the flight. Figure 8 shows the 
mean profdes of 532 nm ATB computed using the standard method (blue) averaged to 
the resolutions of the HSRL products, as well as the ATB using the HSRL method (red) 
for the grey shaded box in Figure 7b centered around 22:32:22 UTC. Both ATB profiles 
follow the modeled molecular profile closely above the cirrus layer and show similar 
structure inside the cirrus layer. The standard ATB retrieval is about 10 percent higher 
than the ATB computed using the HSRL method within the cirrus layer. This difference 
is likely due to the errors in the calibration technique used in both retrievals. The error in 
the CPL calibration constant is estimated to be around 5 percent at 532 nm due to signal 
noise and the presence of aerosols in the CPL calibration zone (Campbell et al., 2008; 
Vaughan et al. 2010). Errors in the determination of the etalon defect parameter can lead 
to errors of as much as 5 percent in the HSRL retrieved attenuated molecular and 
particulate backscatter. Although this comparison provides confidence in the ACATS 
HSRL algorithms, it does not resolve any possible instrument biases. To address this 
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issue, the ACATS standard backscatter and HSRL products are compared to coincident 
CPL cloud and aerosol properties during the WAVE campaign in a companion paper. 

5.0 Summary 

A new multi-channel direct-detection Doppler wind lidar has been developed at 
NASA GSFC for use on the NASA ER-2 called the Airborne Cloud-Aerosol Transport 
System (ACATS). ACATS employs a Fabry-Perot interferometer to provide the spectral 
resolution needed to retrieve the Doppler shift, similar to the ground-based University of 
Michigan MC direct-detection Doppler wind lidar (McGill et al. 1997a). The ACATS 
instrument design includes a seeded laser and circle-to-point converter, as well as a 
heating/cooling loop for stable laser performance during airborne operation. The ACATS 
telescope rotates to four look angles to permit the retrieval of the horizontal wind velocity 
within atmospheric layers. ACATS also advances the technology of a MC direct- 
detection Doppler wind lidar by demonstrating the utility of such an instrument for HSRL 
retrievals of cloud and aerosol properties. 

The nature of a MC direct-detection Doppler wind lidar such as ACATS permits 
three types of cloud and aerosol lidar retrievals: standard backscatter lidar products such 
as ATB and layer boundaries, directly retrieved cloud and aerosol optical properties such 
as extinction and lidar ratio using the HSRL technique, and horizontal wind velocity of 
the cloud or aerosol particles within an atmospheric layer. This paper outlines the 
retrieval algorithms for all two of these types of ACATS data products, focusing on the 
HSRL derived cloud and aerosol properties. The first ACATS science flights were 
conducted during the WAVE project in September 2012. Initial results demonstrate the 
effectiveness of ACATS as an airborne HSRL system. The HSRL ATB retrieval for 
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cirrus observed during the 14 September flight at the 270-degree look angle agrees with 
the ATB derived using the standard backscatter method to within 10 percent. Since the 
ISS CATS HSRL receiver is designed similar to ACATS, the algorithms and data 
products developed for ACATS have direct application to this future spacebome mission 
Furthermore, the ACATS HSRL and wind products can be used for science applications 
such as aerosol transport, smoke plume properties and convective outflow in tropical 
storms. 
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Table 1. Primary system parameters for AC ATS lidar. 

Parameter Value 


Laser Type 
Wavelength 
Laser Repetition Rate 
Laser Output Energy 
Telescope Diameter 
Viewing Angle 
Telescope LOV 
Bandpass Lilter 
Etalon Spacing 
Etalon Reflectivity 
Orders Imaged 
Detector Channels 
Raw Range Resolution 
Horizontal Resolution 
Platform Speed 
Platform Altitude 


Nd: YAG, seeded 
532 nm 
250 Hz 
~10 mJ/pulse 
8 inches 
45 degrees 

350 uradians (full angle) 
150 pm FWHH 
10 cm 
85% 

1.2 

24 

30 m 

1 sec (~200 m) 
~200 m/s 

~ 20 km (65,000 ft) 
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Table 2. Definitions of parameters found in the ACATS lidar equation 


Variable 

Definition 

Units 

N(r) 

number of photons detected per range bin 

- 

r 

distance to the scattering particle 

m 

j 

detector channel 

- 

Et 

transmitted laser energy 

J 

X 

laser wavelength 

m 

h 

Planck's constant 

J sec 

c 

speed of light 

m s' 1 

Aj 

area of lidar telescope 

m 2 

Ar 

range bin width 

m 

Qe 

detector quantum efficiency 

- 

T 0 

system optical efficiency 

- 

T F 

optical filter efficiency 

- 

0 A (r) 

overlap function 

- 

n c 

number of detector channels 

- 

Oc 

detector normalization 

- 

N F sr 

free spectral range (channel number) 

- 

AXpsR 

free spectral range (wavelength) 

m 1 

ax l 

laser broadening 1/e width (wavelength) 

m 1 

AX m 

molecular broadening 1/e width (wavelength) 

m 1 

a(r) 

attenuated particulate backscatter coefficient 

m 1 sr 1 

w(r) 

attenuated molecular backscatter coefficient 

m 1 sr 1 

cp 

off-nadir pointing angle 

degrees 

UlOS 

LOS wind velocity 

m s' 1 

K 

center position of the laser linewidth 

m 

K 

center wavelength of the etalon 

m 
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Figures 


HSRL Total Measured Spectrum ACATS Totol Measured Spectrum 



Figure 1 . The ACATS method images the grey shaded area of the returned atmospheric signal (a) onto a 24 channel 
array detector, which measures the photon counts at each wavelength interval independently as a total backscattered 
signal (b). 
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ACATS Total Measured Spectrum 


ACATS Wind Signal 





Figure 2. The Doppler shifted atmospheric signal (purple) measured by ACATS is compared to an unshifted reference 
spectrum (a), which yields the Doppler wind signal (b) of the ACATS measurement. The broad Rayleigh scattered 
spectrum (c) is measured by ACATS as a nearly flat background of the total atmospheric return signal, resulting in a 
sharp particulate spectrum (d) that is directly measured. 
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652 

653 Figure 3. The fully assembled ACATS instrument (a) includes the receiver tube covered in insulation (left) and a 

654 pressurized telescope dome (right). A picture of the inside of the receiver subsystem (b) shows the etalon (silver device 

655 in the middle), the 24-channel array detector, and circle-to-point converter. The inside of the telescope subsystem (c) 

656 contains a motor to rotate the telescope and a HOE. 

657 
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ACATS Optical Schematic 


si 
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663 



Telescope 

Figure 4. The ACATS optical schematic shows the outgoing 532 nm laser light (dashed green), originating from the 
Nd:YAG laser, directed out of the telescope by a mirror. The return signal (solid green) is passed through the telescope 
and into the receiver subsystem using an optical fiber, where it is transmitted through optical lenses and filters, 
including the etalon. 
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ACATS Corrected Counts for All Look Anales: 26SeDl2 



00:29:20 00:30:28 00:31:38 00:32:45 00:33:55 

Time (UTC) 


665 Figure 5. ACATS photon counts from an ER-2 flight on 26 September 2012. The high count rates between 10 and 12 

666 km show the detection of a cirrus layer at all four look angles at intervals of 60 seconds. 
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a) 


ACATS 532 nm Attenuated Total Backscatter: 26 Sep. 2012 

Resolution: 40 m (vertical), 400 m (horizontal), 000 Degree Look Angle 
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b) ACATS 532 nm Attenuated Particulate Backscatter: 26 Sep. 2012 


Resolution: 450 m (vertical), 5 km (horizontal), 000 Degree Look Angle 
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c) ACATS 532 nm Particulate Extinction Coefficient: 26 Sep. 2012 
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Figure 6. The ACATS 532 nm Attenuated Total Backscatter (km" 1 sr" 1 ) computed using the standard method (a), the 
Attenuated Particulate Backscatter (km" 1 sr" 1 ) derived using the HSRL method (b), and the directly-retrieved Particulate 
Extinction Coefficient (km" 1 ) at the 0 degree look angle (c) for the ER-2 flight on 26 September. ACATS observed 
clouds as the ER-2 flew over the Ohio River Valley (20:28:05 to 21:30:00 UTC) and a large smoke plume (00:24:10 to 
02:10:00 UTC) that extended as high as 6 km over Montana. 
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ACATS 532 nm Attenuated Total Backscatter: 14 Sep. 2012 
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Resolution: 40 m (vertcal), 400 m (horizontal), 270 Degree Look Angle 
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Figure 7. The ACATS 532 nm Attenuated Total Backscatter computed using the standard method (a) and using the 
HSRL method (b) at the 270 degree look angle for the ER-2 flight on 14 September. The grey box focuses on a 35 
minute segment in which the mean profiles are compared in Figure 8 for cirrus clouds. 



ACATS 532 nm HSRL Attenuated Total Backscotter: 14 Sep 2012 
Resolution: 450 m (vertical), 5 km (horizontal), 270 Degree Look Angle 
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ACATS Backscatter (270 deg): 14 Sep. 2012 
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681 Figure 8. The ACATS mean profiles of 532 nm Attenuated Total Backscatter computed using the standard method 

682 (blue) averaged to the resolutions of the HSRL products, as well as the Attenuated Total Backscatter using the HSRL 

683 method (red) for the grey shaded box in Figure 7b (22:1 1:43 - 22:46:21 UTC). 
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